The onset of classical QCD dynamics in relativistic heavy ion collisions 
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The experimental results on hadron production obtained 
recently at RHIC offer a new prospective on the energy de- 
pendence of the nuclear collision dynamics. In particular, it is 
possible that parton saturation - the phenomenon likely pro- 
viding initial conditions for the multi-particle production at 
RHIC energies - may have started to set in central heavy ion 
collisions already around the highest SPS energy. We exam- 
ine this scenario, and make predictions based on high density 
QCD for the forthcoming = 22 GeV run at RHIC. 



High energy nuclear collisions allow to test QCD at 
the high parton density, strong color field frontier, where 
the highly non-linear behavior is expected. Already after 
one year of RHIC operation, a wealth of new experimen- 
tal information on multi-particle production has become 
available It appears that the data on hadron mul- 

tiplicity and its energy, centrality and rapidity depen- 
dence so far are consistent with the approach |^,^ based 
on the ideas of parton saturation |0,|||and semi-classical 
QCD ("the color glass condensate" ) ||9|JlO|] . The central- 
ity dependence of transverse mass spectra appears to be 
consistent with this scenario as well [ |ll| . 

Strictly speaking, the use of classical weak coupling 
methods in QCD can be justified only when the "sat- 
uration scale" Ql [RJy], proportional to the density of 



the partons, 
as{Ql) « 1. 



becomes very large, Qi ^ 



A2 



and 



At RHIC energies, the value of satura- 
tion scale in Au — Au collisions varies in the range of 
Ql = 1 2 GeV^ depending on centrality. At these 
values of Q^, we are still at the borderline of the weak 
coupling regime. However, the apparent success of the 
saturation approach seems to imply that the transition 
to semi-classical QCD dynamics takes place already at 
RHIC energies. 

This may shed new light on the mechanism of hadron 
production at lower energies, perhaps including the en- 
ergy of CERN SPS. Indeed, extrapolating down in en- 
ergy using the formulae of j|] yields for saturation scale 
in central Pb — Pb collisions at SPS energy of = 17 



GeV the value of 1.2 GeV^ |] . The same average 

value at a RHIC energy of ^/s = 130 GeV is reached in 
peripheral Au — Au collisions at impact parameter & w 9 
fm and an average number of participants of Npart ~ 90. 
At Npart < 100, and impact parameters 6 > 9 fm, re- 
construction of the geometry of the collision and the ex- 
traction of the number of participants face sizable un- 
certainties, and no firm conclusion on the applicability 
of the saturation approach can be drawn from the data. 
Given this uncertainty, one may consider two different 
scenarios: 

1) the onset of saturation occurs somewhere in the 
RHIC energy range, below ^/s = 130 GeV but above 
^/s — 17 GeV; the mechanisms of multi-particle produc- 
tion at RHIC and SPS energies are thus totally different; 

2) saturation sets in central heavy ion collisions al- 
ready around the highest SPS energy. The second sce- 
nario would, in particular, have important implications 
for interpretation of the SPS results. 

It should be possible to distinguish between these two 
scenarios by extrapolating the results of Refs. down 
in energy and comparing them to the data. In fact, very 
soon RHIC will collect data at the energy of y/s — 22 
GeV, not far from the highest SPS energy of ^/s = 17 
GeV. In this Letter, we make predictions for hadron pro- 
duction at this energy based on the saturation scenario. 
It should be stressed that a priori there is no solid rea- 
son to expect this approach to work at low energies; we 
provide these predictions to make it possible to decide 
between scenarios 1) and 2) listed above basing on the 
data when they become available. 

In Ref. 1^ we derived a simple analytical scaling for- 
mula, describing the energy, centrality, rapidity, and 
atomic number dependences of hadron multiplicities in 
high energy nuclear collisions: 



^We use Ql oc s^/^ with A « 0.25 0.3 as it follows from 
the scaling behavior of the HERA data [^; see below. 
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with Ql{s) — QI{sq) (s/so)'^/^. Once the energy- 
independent constant c ~ 1 and (^^(so) are determined 
at some energy so, Eq. contains no free param- 

eters. (The value of A, describing the growth of the 
gluon structure functions at small x can be determined 
in deep-inelastic scattering; the HERA data are fitted 
with A ~ 0.25 0.3 ki y ^ the expression @j 

coincides with the one derived in and extends it to 
describe the rapidity and energy dependences. 

Using the value of Ql ~ 2.05 GeV^ extracted at 
= 130 GeV and A = 0.25 |l| used in §, equation (0) 
leads to the following approximate formula for the energy 
dependence of charged multiplicity in central Au — Au 
collisions: 
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At = 130 GeV, we estimate from Eq.(|) 

2/Npart dNch/drj |,,<i= 3.42±0.15, to be compared to the 
average experimental value of 2/Npart dNch/drj |rKi~ 
3.37 ± 0.12 |l|-|. At = 200 GeV, one gets 
2/Npart dNch/drj \n<i= 3.91±0.15, to be compared to the 
PHOBOS value § of 2/Npart dNch/drj |^<i= 3.78±0.25. 
Finally, at = 56 GeV, we find 2/Npart dNch/drj |r;<i = 
2.62±0.15, to be compared to ^ 2/Npart dNch/drj |,,<i = 
2.47 ± 0.25. Having convinced ourselves that our result 
(^ describes the experimentally observed energy depen- 
dence of hadron multiplicity in the entire interval of ex- 
isting measurements at RHIC within error bars, we can 
extrapolate it to the small energy of y/s = 22 GeV and 
make a prediction: 

_^dNch\ ^195 ±0.1; yi = 22GeV. (3) 

Npart drj I 

It is also interesting to note that formula (^, when 
extrapolated to very high energies, predicts for the LHC 
energy a value substantially smaller than found in other 
approaches: 



2 dNch 



Npart drj , 



= 10.8 ±0.5; ^/s = 5500 GeV, (4) 



corresponding only to a factor of 2.8 increase in multi- 
plicity between the RHIC energy of y/s = 200 GeV and 



the LHC energy of ^/s = 5500 GeV (numerical calcula- 
tions show that when normalized to the number of par- 
ticipants, the multiplicity in central Au — Au and Pb— Pb 
systems is almost identical). 

Let us now turn to the centrality dependence. Our 
method has been described in detail before |^,^. We 
first use Glauber approach to reconstruct geometry of 
the collision, and then apply semi-classical QCD to eval- 
uate the multiplicity of produced gluons at a given cen- 
trality and pseudo-rapidity. The Glauber formalism 
(see [^,|l^ for details) allows to evaluate the differential 
cross of inelastic nucleus-nucleus interaction at a given 
(pseudo)rapidity rj: 
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dn 



= j d%V{n-b) (l-Po(&)); 



(5) 



Pa{b) is the probability of no interaction among the nuclei 
at a given impact parameter b: 



Poib) = il-<JNNTAB{b)) 



AB 



(6) 



where (Jnn is the inelastic nucleon-nucleon cross sec- 
tion, and TAB{b) is the nuclear overlap function for the 
collision of nuclei with atomic numbers A and B; we 
have used the three-parameter Woods-Saxon nuclear 
density distributions |1^. For ^/s = 22 GeV we use 
cwAf = 33 ± 1 mb basing on the interpolation of ex- 
isting pp data The correlation function V{n; b) has 
a Gaussian form described in [|l3|,H . The total cross sec- 
tion of inelastic hadronic Au-Au interactions computed 
in our approach at \/s = 22 A GeV is atot = 6.9 ± 0.05 
barn. 

The correspondence between a given centrality cut and 
the mean numbers of nucleon participants and nucleon- 
nucleon collisions can now be established by computing 
the average over the distribution (||), as described in 
At = 22 GeV for Au — Au collisions we find for the 



— 6% centrality cut (iVj 



part} 



332 ±2; (TV, 



coll] 



828±6. 



For 15—25% centrality cut, we have (Npart) — 179±2 and 
(Ncoii) — 367±6, while for the 35 — 45% cut, correspond- 
ing to rather peripheral collisions with the average impact 
parameter of (6) « 9.5 fm, one gets (Npart) = 77± 2 and 
(Ncou) = 120 ± 5. 

We now have the information about the geometry 
of the collision needed to proceed with our calculation 
of centrality dependence of hadron multiplicities in the 
semi-classical QCD approach. However, in applying this 
method at small energies and/or for peripheral collisions, 
we face a fundamental dilemma. In semi-classical ap- 
proach, the multiplicity of the produced gluons is pro- 
portional to l/as(Qs) (this is, of course, the origin of 
the factor ln{Ql/ Aq^,j~i) in our formula (|l|)). Once the 
saturation scale becomes small, the result thus be- 
comes sensitive to the behavior of the strong coupling in 
the infra-red region. Taking a conservative viewpoint. 
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this simply signals that the method ceases to be applica- 
ble. If we accept this, we have to stop and conclude in 
favor of scenario 1) described above. 

However, this is not necessarily correct - there is a solid 
body of evidence from jet physics that QCD coupling 
stays reasonably small, {cts)iR = 0.4-^0.6 in the infrared 
region The "freezing" at small virtualities solution 
for the QCD couphng {as)iR ~ 0.43 has been found by 
Gribov as a consequence of "super-critical" screening 
of color charge by light quark-antiquark pairs. Matching 
QCD onto the chiral theory through scale anomaly leads 
to the coupling frozen in the infrared region, with mag- 
nitude {as)iR « 0.56 "Freezing" solutions for the 
running coupling are repeatedly discussed; see j20| for a 
recent review. It is possible that the presence of relatively 
large scale in hadro-production reflects the properties of 
QCD vacuum |^J||. 

We will thus try to adopt an optimistic point of view 
and assume that the strong coupling indeed "freezes" be- 
low « 0.8 GeV^ at the value of {as)iR ~ 0.5. In fact, 
one may even dare to go further - assuming the validity 
of semi-classical QCD approach at low energies, the cen- 
trality dependence of hadron multiplicity may be used to 
glean information about the behavior of strong coupling 
in the infra-red region. 
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FIG. 1. Centrality dependence of the charged multiplicity 
per participant pair at different pseudo-rapidity intervals at 
^/s = 22 A GeV; see text for details. 



To evaluate the resulting centrality dependence around 
ry = we use two different ansatze for the running cou- 
pling: a) "smooth freezing" as ~ 1/ ln((Q^-|-A^)/AQp^), 
with A = 0.8 GeV^; b) "sudden freezing" , when as is sim- 
ply put equal to as(A^) when < A^. The results are 
shown in Fig. ^ by solid (ansatz a)) and dashed (ansatz 



b)) lines. Note that even in the case of "sudden freez- 
ing", centrality dependence is smooth - this is because 
the fraction of the transverse area where the local value 
of becomes smaller than the cutoff A^ is a smooth 
function of centrality. 

Let us now discuss rapidity dependence. Unfortu- 
nately, we have found that at at low energies ^/s ^ 20 
GeV the expression (|l|) provides a poor approximation to 
the numerical result based on the general formula [0jl^] 
used in Q: 



E 



da 



1 

T p2 



dki 



(7) 



where xi^2 — {Pt/Vs) exp(±?/) and ipA{x, k^) is the unin- 
tegrated gluon distribution. This happens because at low 
energies the limited phase space suppresses the transverse 
momentum distribution of the produced gluons already 
below the saturation momentum Qs- We thus have to 
evaluate the integral in numerically (see for the 
list of formulae needed for this computation). 
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FIG. 2. Pseudo-rapidity distribution of the charged mul- 
tiplicity per participant pair in different centrality cuts at 
= 22 A GeV. 

To convert the computed rapidity distributions of 
gluons to the observed pseudo-rapidity distribution of 
hadrons, we follow the procedure of [ p[, assuming that 
the "local parton-hadron duality" (see and references 
therein) in the space of emission angles 9, or, equivalently, 
in pseudo-rapidity rj = — lntan(0/2). This corresponds 
to the physical assumption that once a gluon has been 
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emitted along a certain direction, its final state interac- 
tions and fragmentation will not significantly change the 
direction of the resulting hadrons. The results of our 
calculations are presented in Figs. ^ and 

How reliable are our results, apart from the obvious 
leap of faith involved in the application of semi-classical 
method at small energy? The least reliable of our pre- 
dictions is the distribution in pseudo-rapidity; indeed, 
we have assumed that multi-particle production is domi- 
nated by gluons, and this is not necessarily so at ^/s = 22 
GeV, where valence quarks may give essential contribu- 
tion even at central rapidity. Also, since we are quite 
close to the fragmentation region at this energy, multi- 
parton correlation effects can also play a role. The ab- 
solute value of multiplicity around 77 = is more stable, 
but still may be affected by the contribution from quarks 
and deviations from behavior of the nuclear gluon 

distribution at larger x. 

There does exist however a prediction of our approach 
which is both robust and distinct: it is the rise with cen- 
trality of multiplicity per participant shown in Fig. |l| 
around 77 « 0. The shape of this dependence simply re- 
flects the running of strong coupling, and is similar to the 
shape observed at much higher energy of ^/s = 130 GeV: 
indeed, the logarithm in Eq.(|^) is a slowly varying func- 
tion, and the dependence of saturation scale on energy is 
quite weak. This prediction is in marked contrast both to 
the strong increase of the slope of centrality dependence 
with energy predicted in a two-component model p^ ] (for 
a recent development, see however [Q) and to the final- 
state saturation model p6[ |, predicting a nearly constant, 
weakly decreasing, centrality dependence of multiplicity 
per participant. 
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